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Abstract 

General asymptotic formulas are given for the coefficient Cg of the term of 
multipole number i in the temperature correlation function of the cosmic mi- 
crowave background, in terms of scalar and dipole form factors introduced in 
a companion paper. The formulas apply in two overlapping limits: for £ 3> 1 
and for id/dA <^ 1 (where dA is the angular diameter distance of the surface 
of last scattering, and c? is a length, of the order of the acoustic horizon at 
the time of last scattering, that characterizes acoustic oscillations before this 
time.) The frequently used approximation that Cg receives its main con- 
tribution from wave numbers of order i/dA is found to be less accurate for 
the contribution of the Doppler effect than for the Sachs-Wolfe effect and 
intrinsic temperature fluctuations. For id/dA -C 1 and i >2, the growth of 
Ci with i is shown to be affected by acoustic oscillation wave numbers of all 
scales. The asymptotic formulas are applied to a model of acoustic oscilla- 
tions before the time of last scattering, with results in reasonable agreement 
with more elaborate computer calculations. 
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I. INTRODUCTION 



A companion paper[l] has shown how to express the temperature fluc- 
tuation in the cosmic microwave background in any direction as an integral 
involving scalar and dipole form factors F{k) and G{k), which characterize 
acoustic oscillations before the time of last scattering. In the present paper 
we derive asymptotic formulas for the strength Ce of fluctuations at multi- 
pole number £ for form factors of arbitary functional form. After outlining 
our assumptions and reviewing some generalities in Section II, our general 
result in the limit of £ ^ 1 [Eq. (26)] is derived in Section III. In this limit 
£{£+ l)Ci depends on £ and the angular diameter distance dA at the time of 
last scattering only through the ratio i/d-A- (This is why the heights of the 
Doppler peaks do not depend on parameters like the cosmological constant 
that affect dA but not the form factors.) Our result in the limit £d/dA <S 1 
[Eq. (43)] is derived in Section IV. (Here d is some length characterizing 
acoustic oscillations, such as the acoustic horizon distance dn at the time of 
last scattering) . These ranges of £ overlap because dA ^ d. 

Even without a detailed calculation of the form factors, these results have 
a moral for the physical interpretation of measurements of C^. It is common 
to interpret these measurements by supposing that Ce arises mostly from 
fluctuations of wave number k ~ £/dA- Eq. (27) shows that this is a fair 
approximation for the contribution of the scalar form factor F{k), which 
represents the Sachs-Wolfe effect and intrinsic temperature fluctuations; Q 
receives no contribution from F{k) with k < £/dA, and the contribution from 
k 3> i/dA is suppressed by a factor /3^^(/3^ — 1)~^/^, where /? = kdA/ £■ In 
particular, a peak in the magnitude of the scalar form factor F{k) at some 
wave number ki (like the peak found in the simple model studied in reference 
[1] at k — n / dn) will show up in £(£-|- 1)Q at a value of £ less than but close 
to kidA- For instance, we will see in Section V that the peak in |-F(A;)| at 
k = ir/dn produces a peak in £{£ + 1)C^ at £ ~ 2.6dA/dH rather than at 
TTdA/dn- But Eq. (27) also shows that this interpretation of Q is much less 
useful for the contribution of the vector form factor G{k), which arises from 
the Doppler effect; Ce also receives no contribution from G{k) with k < £/dA, 
but instead of the contribution from k ~ £/dA being enhanced by a factor 
(/9^ — 1)~^/^, it is suppressed by a factor (/3^ — 1)-*^/^. Indeed, we will see 
in Section V that for sufficiently small baryon number the peak in G{k) at 
k = Tx/2dH found in the simple model of reference [1] does show up as a peak 



1 



in £{£+l)C£, but at £ ~ AbdA/dn, much less than [tt /2)dA/ du- Furthermore, 
the behavior of + 1)Q for IdjdA near zero depends on the values of 
and G{}z) for all k. This points up the value of observations that can measure 
the correlation function of temperature fluctuations directly, as a supplement 
to measurements of C^. 

The results obtained in Sections III and IV are used in Section V to 
calculate Ci for the approximate form factors calculated in reference 1. In 
agreement with what is found in more accurate computer calculations, the 
position l\ of the first Doppler peak is not a sensitive function of the baryon 
density parameter O^/i^. On the other hand, we find that the ratio of the 
value of i{JL + 1)C^ at the first Doppler peak to its value at £ <C dA/dn is a 
sensitive indicator of the value oiVLsh^- 



II. GENERALITIES 



The companion paper[l] shows that, in very general models (but assum- 
ing only compressional normal modes, with no gravitational radiation), the 
fractional variation from the mean of the cosmic microwave background tem- 
perature observed in a direction n takes the general form 



Ar(n) 



/ 



d^k ek e 



F{k) + in-kG{k) 



(1) 



aside from effects arising from late times, which chiefly affect the coefficients 
Q for relatively small i. Here dA is the angular diameter distance of the 
surface of last scattering 



1 + zl) 



sinh 



1/2 

C 



dx 



(2) 



where Vie = I — VL^ — ^m^ and Vt^ and VLm are the present ratios of the energy 
densities of the vacuum and matter to the critical density 3Hq/87iG. (If the 
vacuum energy were to change with time, as in theories of quintessence, then 
the formula for dA would need modiflcation, but there would be essentially 
no change in the other ingredients in Eq. (1), as long as the quintessence 
energy density makes a negligible contribution to the total energy density 
at and before the time of last scattering.) Also, k^ek is proportional (with 
a k-independent proportionality coefficient) to the Fourier transform of the 
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fractional perturbation in the energy density early in the radiation-dominated 
era. The averageQ of the product of two es is assumed to satisfy the condi- 
tions of statistical homogeneity and isotropy: 

{e^e^,) = 6'\k + k')V{k) (3) 

with k = |k|. The power spectral function V{k) is real and positive. Where 
a specific expression for V{k) is needed, we will use the 'scale-invariant' (or 
n = 1) Harrison-Zel'dovich form suggested by theories of new inflation: 

V{k) = Bk-\ (4) 

with B a constant that must be taken from observations of the cosmic mi- 
crowave background or condensed object mass distributions, or from detailed 
theories of inflation. 

The form factors F{k) and G{k) characterize acoustic oscillations, with 
F{k) arising from the Sachs-Wolfe effect and intrinsic temperature fluctu- 
ations, and G{k) arising from the Doppler effect. For instance, they are 
calculated in reference 1 in the approximation that perturbations in the 
gravitational field at and before the time of last scattering arise entirely 
from perturbations in the density of cold dark matter. For very small wave 
numbers the form factors are 

F{k) ^ 1 - ?,khi/2 - 3[-r' + in(i + + • • • , (s) 

G{k) 3ktL - 3Pti/2(l + + • • • , (6) 



while for wave numbers large enough to allow the use of the WKB approxi- 

ktL > e (7) 



mation, i. e.. 



the form factors are 

F{k) = (1 - 2^/kHl)-^ [-3^ + 2^/kHl + (1 + ^)-i/4g-fc2rfz, cosikdn)] , (8) 
and 

G{k) = 73 (1 - 2^/kHl)-\l + ^)-3/4e-fc^rfi, smikdn) . (9) 

**The average here is over an enserable of possible fluctuations. Using Eq. (3) to analyze 
the particular element of this sample observed in our universe relies on ergodic arguments, 
which are not exact except in the limit £ — s- oo. However, corrections are manageable[2] 
even for small £. 
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Here ti is the time of last scattering; ^ is 3/4 the ratio of the baryon to 
photon energy densities at this time: 



'3p 



B 



(10) 



dn is the acoustic horizon size at this time, and d^) is a damping length, 
given by Eq. (48). These formulas for the form factors are mentioned at 
this point only for illustration; wc will be working here with general form 
factors F{k) and G{k), and will not make use of the specific formulas (5)- 
(10) until Section V. But we will assume throughout that any lengths d that 
(like du and do in Eqs. (8) and (9)) characterize the /c-dependence of the 
form factors are much smaller than the angular diameter distance of last 
scattering. This is a good approximation: for instance, if the ratios of matter 
and vacuum energy densities to the critical density have the present values 
Qm = 0.3 and Qa = 0.7, then dA/dn runs from 91.7 to 79.7 for values of 
VLsh? running from zero to 0.03, and do is smaller than dui independent of 
the value of Hq. 

It is usual to employ the well-known expansion of a plane wave in Legen- 
dre polynomials, and write Eq. (1) as 



Ar(n) 
T 



+j'e(kdA) G(k) 



n 



k) 



3e{kdA) F{k) 



Using Eq. (3) and the orthogonality property of Legendre polynomials 

47r 



j dnj^Pe(n-k) Pe (n'-k) 
one finds that 



2£+l 



5a'Pi{n-h') , 



'Ar(n) AT(n')' 



T 



T 



ff^l C,P,{n-n') , 



£=0 



with the conventional coefficient Ci taking the value 



Q = 16r 



/■oo 

Jo 



k^ dk 



jeikdA) F{k) + j[{kdA) G{k) 



l2 



(11) 



(12) 



(13) 



(14) 
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This familiar formula is adequate for numerical calculation of Ce, but it hides 
the essential qualitative aspects of the dependence of Q on £: that Ce for 
£ ^ 1 depends on the ratio i/dA, and that £{i + 1)Q approaches a constant 
for sufficiently small values of this ratio, whether i itself is large or small. To 
obtain these results, we must now distinguish between the two cases i ^ 1 
and £ ^ rf^/rf (but £ > 2), where d is a typical length characteristic of the 
form-factors F{k) and G{k). These two cases overlap because, as remarked 
above, dA is much larger than d. 

III. LARGE £ 




The usual way of obtaining the contribution of the scalar form factor to 
Ce for large £ is to note that the integral (14) receives its largest contribution 

when the argument of the spherical Bcssel function is of order £, in which 
case we can use the approximation that, for £ ^ oo, 

z < u 

- i/^)'~-'^/^cos (^Vz"^ — i/"^ — 1/ Bxccos{v / z) — \) z > v , 

(15) 

where z/u is held fixed at a value ^ 1, with u = £ -{- 1/2. The procedure is 
straightorward for the terms in Eq. (14), but for the FG and terms 
involving the Doppler effect we run into a difficulty: differentiating the factor 
(z^ — i/^)~^/^ in Eq. (15) yields larger negative powers of z^ — v^ that introduce 
divergences from the part of the integral in Eq. (14) near the lower bound 
k = v/dA- These infrared divergences are spurious, because the asymptotic 
formula (15) breaks down if we let z and v go to infinity in such a way 
that z/y — > 1. This problem can be dealt with by switching to a different 
asymptotic limit [3] for k near v/dA- Here we will use a different method [4] 
which avoids the delicate problem of the asymptotic behavior of j(.{z) and 
j[{z) for z near v. 

We return to Eq. (1), and use Eq. (3) to put the correlation function of 
observed temperature fluctuations in the form 



'AT(n) AT(n') 



j d^kV{k) exp (idAk ■ (n - n')) [F'^{k) 



T T I 

+ik-{n-h')F{k)G{k) + {k-h){k-h')G^{k)] . (16) 
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The integral over the direction of k is easy, and gives the correlation function 



/AT{n) AT(n') 
1 



T 



POO 

Att / k^dkV(k) 

Jo 



F\k) + F{k)G{k) 



d 



d{dAk) 



sm{dAk9) 
dAkO 



(17) 



where 9 = \n — n'\. (This formula may prove useful in analyzing observations 
that give the correlation function directly, rather than in terms of C^.) The 
amplitude Q is defined as the integral 



(18) 



where ^ = h ■ h' = 1 — 9"^ /2. For large i the Legendre polynomial P^(yu) 
oscillates rapidly for 9 ^ so the integral is dominated by values of 9 
of order in which case we can use the well-known limiting expression 
Pi{lj) Jo{^9), and write 



Ce SttM k^dkV(k) JqU9)9 
Jo Jo 



d9 



F'^{k)+F{k)G{k) 
39^\ 



d 



did^kf 



d{dAk) 
sai{dAk9) 



dAk9 



■ (19) 



The integral over k is dominated by values for which kdA9 is of order unity, 
so the derivative d/d{dAk) is effectively of order 9 ^ l/£. Thus to leading 
order in l/£, Eq. (19) may be simphfied to 



POO p2 

Ce SttM k^dkr{k) Jo{i9)9 
Jo Jo 
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d9 



sm{dAk9) 
dAk9 



(20) 



Introducing a new variable s = £9 and changing the upper limit on the 
s-integral from 2£ to infinity, we may write this as 

S7]-2 POO poo 

Ce ^ —j^ k^dkV{k)J^ Ms)sds 

sm{dAks/i) 



X 



(dAks/i) 



(21) 



The integral over s is easy for the term; we need only use the formula[5]: 



Jo(s) sin(/?s) ds 



{(3^ - 1)-V2 



/? < 1 
P>1, 



(22) 



where here /3 = dAk/i. The integral of the term takes a httle more work. 
We use the formula (1 + d? /dx^) svax/x — —{2/x)d/dx{smx/x) and do the 
remaining integral by parts, so that 



/ Ms) 

Jo 



1 + 



92 



sin(/3s) , 2 



(</2('S) + -^o("S)) sin(/3s) 
Here we also need the formula[5] 



roo 

^0 



d sin(/3s) 
9s s 



(is 



2 1 
^ ~ ^3 Jo 



(23) 



/■oo 

/ J2(s) sin(/9s) ds 

JO 

so that 

Jo(s) s 



2/3 

-(/32-i)-V2(/3 + VF^)-^ 



/5 < 1 
/3 > 1 ' 



1 + 



92 



5(/3s)^ 



sin(/3s) 



ds 







P < 1 
/3 > 1 . 



(24) 



(25) 



Using Eqs. (22) and (25) in Eq. (21) then gives our final general formula for 
Q at large i: 



Sir^i /-oo 



d\ 



(3F\ip/dA) Vp^^G\i(3/dA) 



+ 



■ (26) 



Note that i'^Ci depends on i and d^ only through its dependence on the ratio 
i/dA. 

For instance, if take the power spectral function to have the scale- invariant 
form V{k) = Bk-^, then for £ > 1 



dp 



F\ip/dA) , VP^^G\ep/dA) 



+ 



P' 



(27) 



(We have taken advantage of the fact that here we are considering £ ^ 1 
to change a factor £^ to £(£ +1), in order to facilitate comparison with the 
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results of the next section.) The rapid fall-off of the coefficient of F'^ for 
/5 > 1 suggests that the contribution of the scalar form factor F to Q is 
dominated by wave numbers close to d^/^, as is usually assumed. On the 
other hand, the contribution of the dipole form factor G{k) for wave numbers 
immediately above is actually suppressed by the factor \/ (5'^ — 1 in the 
second term of Eq. (27). 

IV. SMALL ^d/dA 

Here we will adopt the 'n = 1' scale-invariant spectrum Vik) ~ Bk~^ 
from the beginning, so that the general formula Eq. (14) becomes 



mFi-\+j',{s)G{- 



di 



ds 



(28) 



To generate a series for £(£ -|- 1)C^ in powers of £/o?a we expand the form 
factors in power series: 



F{k) = Fo + F2 A;' + 



G{k) = Gik + G^k^ + 



(29) 



(The power series for F and G must be respectively even and odd in /c, in 
order that the integrand in the temperature fluctuation (1) should be analytic 
in the three- vector k at k = 0.) The leading term in is well known; using 
a standard formula[6] : 



/■oo 

/ m 

Jo 



2^-\r{2 - m) r -t- 



P2 ^3^^ r + 2 - 

we find the term in Eq. (28) of zeroth order in 1/dA- 



(30) 



Sn^BFl 

e{i + i) 



(31) 



There is no difficulty in also calculating the term in Eq. (28) of first order in 

l/dA-. 



^(1) ^ / 327r^^FoG'/ 



je{s)ji{s) ds 



IGtt^ BFqGi 
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00 
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(32) 
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But we run into trouble in calculating the term of second order in 1/dA- The 
second derivative of Cn with respect to 1/dA is 



d{lldAf 



16n'B {j]{s)F'"{s/dA)+J%)G'"{s/dA) 



+2Us)f,is) F{s/dA)G{s/dA) 



s ds 



(33) 



The term doesn't contribute to the part of of second order in 1/dAi 
because F{k)G{k) contains only odd powers of k. To calculate the contri- 
bution of the i'1 term, we need to supplement Eq. (30) with the additional 
formula:[| 



j7(s) s'^-^ds 



1 + 



(m-3)(m-2)((m-2)(m-3) - 2£(£ + 1) 



2(3 



m] 



+ ^ - 1 

^2 



f + 2 



(34) 



The second derivative (33) is divergent at l/ciA = 0, as shown by the factors 
r(2 — m) in Eqs. (30) and (34), which become infinite for m = 2. Of course, 
there is no infinity in Cf, it is simply not analytic in 1/dA at l/rf^ = 0. 

We can deal with this problem by a method similar to the dimensional 
regularization technique used in quantum field theory[7]. We treat m as a 
complex variable that approaches m = 2. In this limit, Eqs. (30) and (34) 



give 



oo 1 
I 



m 



+ J2--C + \n2-D 



r=l 



j'/{s)s"'-^ds 



m — 2 



+ + \n2~ D + 1 



r=l 



(35) 



(36) 



where C is the Euler constant C = -r'(l) = 0.57722, and D = -r'(l/2)/r(l/2) 
1.96351. The important point here is that the parts of the integrals (35) and 
(36) that are divergent at m = 2 are independent of i, and so also is the part 

t This formula was obtained by using the Bessel differential equation to show that 
j'/{z) = z'^)jj{z) + {zjj{z))" I2z, and then using Eq. (30) with two integrations 

by parts. 
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of Ce that is non-analytic in 1/dA at 1/dA = 0. Using Eqs. (29), (35) and 
(36) in Eq. (33) thus gives the part of Ce that is of second order in 1/dA as 

e 1 

CP ^ -Sn^BdA^ (2FoF2 + Gl)j2- + £ -independent terms . (37) 



r=l 



We can check the consistency of these results and calculate the ^-independent 
terms here by using our previous result (27) in the case where i is large and 
id/dA is small, where d is whatever length characterizes the /c-dependence of 
the form factors. The term in Eq. (27) of zeroth order in £d/dA is 



in agreement with Eq. (31). Also, Eq. (27) has no terms of first order in 1/dA, 
in agreement with Eq. (32). To calculate the terms in Eq. (27) of second order 
in 1/dA, we express F'^{k) and G'^{k) in terms of cosine transforms 

F'^{k) = FqV da f(a) (l-cos(ka)) , G^(k) = da g{a) (l-cos(A;a)) . 

(39) 

Then for £ > 1 and id/dA < 1, Eq. (27) gives 



{2F0F2 + GI) (^\n(^- 



id 



+ C 



+ Gl 



(40) 



2dAl ' 2 

where J is a typical value of the variable a in the cosine transforms (39): 

, T_ /o°° [/(«)+ Inada 

Io^[f{a)+g{a)]a^da ^ ^ 

Eq. (40) agrees with the limit of Eq. (37) for large i, because in this limit 
Y,i l/r \ni + C, and now fixes the ^-independent terms in Eq. (37) so 
that, for any i with id/dA <C 1, 

{2FoF, + Gl) (ln( — I+E--0 



G. 



(2) 



d\ 



2d, 



r=l 



(42) 



Putting together Eqs. (31), (32), and (42) gives our final formula for Q 
in the case id/dA ^ 1 and i >2: 



+ l)Gi = Sn^BF^ 1 1 



e{i+i) 

d\ 



d^ In 



Ma. 



d 



+ . . . 



(43) 
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where now we introduce a pair of characteristic lengths: 

2F0F2 + Gl 



r/2 



3F0F2 + jGf 



(44) 



The logarithm in Eq. (43) is large and negative, so i{i + 1)C^ will increase 
or decrease with £ for sufficiently small i according as c?^ > or d"^ < 0. 
(Taken literally, Eq. (43) would suggest that this behavior is reversed when 
the sum over r becomes large enough to cancel the logarithm, but this is at 
£ ~ 2e~'^dA/d^ which is large enough to invalidate the approximations that 
led to Eq. (43).) Note that, while d and d' depend only on the behavior of the 
form factors near zero wave number, the length d given by Eq. (41) depends 
on the behavior of the form factors at all wave numbers. Consequently, 
although the value of at low I depends only on the form factors at /c = 0, 
somewhat surprisingly the growth of Q for small £ depends on the form 
factors at all wave numbers. 

V. APPLICATION 

To illustrate the use of the asymptotic formulas obtained here, we will 
now apply them to the simplified model described in reference 1: the uni- 
verse before last scattering consisting of pressureless cold dark matter and 
a photon-nucleon-electron plasma; no gravitational radiation; and negligible 
contributions of the plasma and neutrinos to the gravitational field. In this 
case, the comparison of Eqs. (5) and (6) for the long wavelength limit of the 
form factors with Eq. (29) gives 



(45) 



so the lengths (44) are here 

d'^ = 6tl , d'^ = . (46) 
Hence Eq. (43) then gives the behaviour of Q for £d/dA <S 1 and £ >2 as 



£{£+l)Ce^87r^B{l- 



U{£ + \)t\ 

"A 



In 



2d, 



+ E 

r=l 



+ 



(47) 



Aside from its weak dependence on rf, the behaviour of Ci for Id/dA ^ 
1 is independent of the baryon density, in agreement with more accurate 
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computer calculations [8]. We can't calculate the length d without a model 
that would give the form factors at all wave numbers, but d is expected to be 
roughly of order dn-, and since dA/dn is large the logarithm is not sensitive 
to the precise value of d. If for instance we take d = \/3tL = dA/58.5 (the 
acoustic horizon at last scattering for Qm = 0.4, Qy = 0.6, and Qb = 0) then 
the quantity £{£+ l)Ce/8TT^B rises from unity when extrapolated to £ = to 
1.044 at £ = 5, and to 1.118 at £ = 10, which is probably the highest value 
of i for which the approximations leading to Eq. (47) are rehable. 

For i of the order of dA/dn the coefficients can be calculated under 
the simplifying assumptions of this section by using the form factors given 
by Eqs. (8) and (9) in Eq. (27). The damping length is given in reference 1 
as 

dl = Vl + AVl ^ 0.029 (y^(^ + 2{i%y ) + ■ ^^^^ 

Our results for Ci at and below the first Doppler peak are not sensitive to 
do- We will simplify our calculations here by dropping the terms in Eqs. (8) 
and (9) that are proportional to the ratio i/kH\, on the grounds that these 
terms are not very different from corrections to the WKB approximation that 
are not included either. (At the first Doppler peak ^/k'^t\ increases with ^ 
and hence with ilsh'^, and for Qsh^ = 0.03 it has the value 0.20. But to 
be honest, the real reason for dropping these terms is that they spoil the 
agreement of our results for the height of the first Doppler peak with more 
accurate numerical calculations.) The results obtained now depend critically 
on the baryon density parameter ~ 27 0^^^, and are shown in Figure 1 for 
values of Qsh^ ranging from zero to 0.03. 

For fifi = (in which case the WKB approximation is not needed, so 
that Eq. (27) should give Ce down to values of i of order two) the behavior 
of Ci is nothing hke what is observed: i{i + l)Q/87r^S rises from unity to 
1.1 at a 'zcroth Doppler peak' at Idu/dA — 0.45 (due to the maximum in the 
Doppler form factor G{k) at kdu = /t/2), then dips to 0.7 at idn/dA — 1.6, 
and then rises again to a first Doppler peak at idn/dA — 2.83. 

For Qfi/i^ > 0.01 the behavior of Ce within the range of vahdity of the 
WKB approximation is much more like what is observed: i{£ + 1)Q rises 
monotonically to a first Doppler peak at idn/dA very roughly of order tt 
(though actually around 2.6). There is another clear peak at £ ~ 8.7 dn/dA, 
presumably arising from the peak in F{k) at k = 3TT/dH- The weaker peaks 
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Figure 1: Plots of the ratio of the multipole strength parameter ^{^+ 1)C£ to 
its value at small versus idn/dA, where dn is the horizon size at the time of 
last scattering and dA is the angular diameter distance of the surface of last 
scattering. The curves are for Qsh^ ranging (from top to bottom) over the 
values 0.03, 0.02, 0.01, and 0, corresponding to C, taking the values 0.81, 0.54, 
0.27, and 0. The solid curves are calculated using the WKB approximation; 
dashed lines indicate an extrapolation to the known value at small idn/dA- 
These results are independent of the parameters Hq, Q\, and Qm- 
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in i{i+l)Ce arising from peaks in F{k) near even values of kdn/TT are absent 
here, presumably because of our neglect of the contribution of radiation and 
neutrinos to the gravitational field. Another difference between the curves of 
Figure 1 and more accurate computer calculations is that, again because we 
neglect the contribution of radiation and neutrinos to the gravitational field, 
our results do not show the fall-off of i{i + l)Ce at large i associated with 
the fall-off of the familiar transfer function T[k) at large k. 

The values of the position iidn/dA of the first peak and the ratio of its 
height £i{ii + l)Cfi to the value Stt^S ~ 6C2 for small i are given for various 
baryon densities in Table 1. These results are independent of other param- 
eters. In the last two columns of Table 1 we also give values of dA/du for 
VLm = 0.3 and fl\ = 0.7, and the corresponding results for the multipole num- 
ber £1 of the first Doppler peak. In calculating the horizon at last scattering 

we have now (somewhat inconsistently) taken into account the effect of 
photons and three flavors of neutrinos and antineutrinos on the expansion 
rate, which gives 



where A = 0.047/fijw/i^ is the ratio of photon and neutrino energy density 
to dark matter energy density at the time of last scattering, and d^ is given 
by Eq. (2). In calculating the values of dA/dn in the table we have taken 
Vtuh? = 0.15. 

We see from Table 1 that the position of the first Doppler peak does not 
depend strongly on Vtsh"^, while its height is a sensitive function of fi^/i^. 
For VLsh^ between 0.02 and 0.03 the height and position are in fair agreement 
with what is observed, though of course the serious comparison of theory with 
observation relies on more accurate computer calculations. The qualitative 
results obtained here suggest that if one were to rely on a single feature of 
the plot of i{l + V)Ci versus £ to measure flsh^, then the ratio of the the 
height of the first Doppler peak to the value for lower i values studied by 
the COBE satellite would be more useful than the ratio of the heights of the 
first and second Doppler peaks, which relies on less precise data, depends on 




(49) 
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Table 1: Location iid^/dfj of the first Doppler peak and height of the peak 
in i{i + l)Ci relative to its value Stt'^B ^ 6C2 for £ extrapolated to zero 
for various values of the baryon density parameter. These results, and the 
curves in Figure 1, are independent of the values of Hq, Q\, and (within our 
approximations) Qm and. The last two columns give the values of dA_/dH 
and ii for Qm = 0.3, = 0.7, with dn calculated taking into account 
the contribution of photons and neutrinos to the expansion rate, and using 
^Mh'^ = 0.15. 







iidn/dA 




dA/du 


£1 








2.83 


0.863 


91.7 


260 


0.01 


0.27 


2.65 


2.34 


87.1 


231 


0.02 


0.54 


2.60 


5.09 


83.6 


217 


0.03 


0.81 


2.58 


9.115 


79.7 


206 



complicated damping effects, and is more sensitive to other parameters, such 
as fijv/^^ and the rate of change, if any, of the vacuum energy. Of course, for 
high precision one must use the whole plot of ^{^ + l)Ci versus C. to measure 
all these parameters together. 
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